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ABSTRACT 

We present a comprehensive grid of synthetic stellar-atmosphere spectra, suitable for 
the analysis of high resolution spectra of hydrogen-deficient post- Asymptotic Giant 
Branch (post-AGB) objects hotter than 50000 K, migrating along the constant lu- 
minosity branch of the Hcrtzsprung-Russell diagram (HRD). The grid was calculated 
with CMFGEN, a state-of-the-art stellar atmosphere code that properly treats the 
stellar winds, accounting for expanding atmospheres in non-LTE, line blanketing, soft 
X-rays, and wind clumping. We include many ionic species that have been previously 
neglected. Our uniform set of models fills a niche in an important parameter regime, 
i.e., high effective temperatures, high surface gravities, and a range of mass-loss values. 
The grid constitutes a general tool to facilitate determination of the stellar parameters 
and line identifications and to interpret morphological changes of the stellar spectrum 
as stars evolve through the central star of planetary nebula (CSPN) phase. We show 
the effect of major physical parameters on spectral lines in the far-UV, UV, and op- 
tical regimes. We analyse UV and far-UV spectra of the central star of NGC 6905 
using the grid to constrain its physical parameters, and proceed to further explore 
other parameters not taken in consideration in the grid. This application shows that 
the grid can be used to constrain the main photospheric and wind parameters, as a 
first step towards a detailed analysis. The full grid of synthetic spectra, comprising 
far-UV, UV, optical, and IR spectral regions, is available on-line. 

Key words: stars: post-AGB - stars: atmospheres - stars: low-mass - stars: mass-loss 
- stars: winds - stars: individual (NGC 6905) 



1 INTRODUCTION 

Evolutionary models predict that low and intermediate mass 
stars (of initial masses between 1 and 8 Mq) leave the AGB 
and evolve through the constant luminosity track of the 
HRD up to very high effective temperatures (~ 200000 
K, for CSPNe of about 0.6 M ), then turn on to the 
White Dwarf (WD) cooling track, where their luminosi- 
ties fade while their effective temperatures decrease and, 



finally, they end their lives as WDs (see, for example, Vassil- 
|iadis and~W ood 1994; Hcrwig 2005[ and references therein). 
As the pos-AGB stars evolve toward higher temperatures, 
they ionize the surrounding hydrogen-rich material previ- 
ously ejected by the AGB precursors, giving rise to bright 
planetary nebulae (PNe). Because these objects expel the 
majority of their initial mass prior to settling on the WD 
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cooling sequence, they are a prime source of chemically pro- 
cessed material for the interstellar medium (ISM) and thus, 
a fundamental ingredient of galactic chemical evolution (e.g. 
|Marigo|200T| |Bianchi et al.|2011a| |Karakas|20T0] and refer- 
ences therein). In addition, their intense UV radiation fields 
and fast winds during the post-AGB phase can influence the 
dynamics of the ISM. 

Approximately 20 per cent of the CSPNe are hydrogen 
deficient, a condition thought to be due to thermal pulses 
that occur after the star leaves the AGB. These pulses can 
happen at different epochs of the evolution, resulting in stars 
with different properties. If the thermal pulse occurs after 
the star has already entered the WD cooling track, it is 
termed a very late thermal pulse or VLTP (from the 'born- 
again' scenario of Iben et al~| l983 ) . The late thermal pulse 
or LTP occurs during the CSPN phase while hydrogen burn- 
ing is still ongoing. In both cases, the star returns towards 
the AGB. The star can also undergo an AGB final ther- 



2 Keller et al. 



mal pulse or AFTP, which happens at the end of the AGB 
phase ( |Blocker|200l| |Herwig||2001[ ). The amount of remain- 
ing hydrogen varies among these scenarios. A VLTP pro- 
duces hydrogen- free stars ( |Werner and Herwig 2006), while 
in the LTP case, the hydrogen content at the star surface is 
decreased to a few per cent by mass. If the star undergoes 
an AFTP instead, it is left with a relatively high hydro- 
gen content. A measurement of the nitrogen abundance can 
help distinguish whether a VLTP or a LTP have occurred, 
because the evolution of the stellar nitrogen content is differ- 
ent in each of these scenarios, resulting in an abundance of 
about 0.1 per cent by mass in the case of a LTP event and up 
to a few per cent in the case of a VLTP (Werner et al.|2006 1. 
The three scenarios also predict somewhat different helium, 
carbon and oxygen abundances, therefore, measurements of 
these quantities help to determine if the star underwent a 



GALEX surveys (Martin et al. 2005 Bianchi 20091, thou 



LTP, VLTP, or a AFTP (see, for example, |Blocker|200l| and 
references therein). 

The H-deficient CSPNe are commonly divided into 
three main classes, which are thought to constitute an evo- 
lutionary sequence. These are: [WC], showing spectra very 
similar to those of Population I Wolf-Rayet (WR) stars, 
with strong carbon and helium emission lines; the PG1159- 
type, occupying the region at the top of the WD cooling 
track in the HRD, and characterized by absorption lines of 
highly ionized He, C and O, besides UV wind lines much 
weaker than the ones seen in [WC] stars; [WC]-PG1159, 
that are believed to be transition objects between the two 
other classes. The [WC] stars are further divided into early 
([WCE]) and late-type ([WCL]) objects, showing lines from 
ions of higher and lower ionization stages, respectively. The 
evolution of the H-deficient CSPNe might proceed in the 
following way: AGB -» [WCL] -> [WCE] — » [WC]-PG1159 
->■ PG1159 -> WD ([Werner and Heber|p9l] |G6rny and 
Tylenda 2000 , Pena et al.|2001| ). According to this scenario, 
the [WC] CSPNe would evolve from the AGB at an almost 
constant luminosity, towards higher temperatures. As the 
stars evolve, their radii decrease until the nuclear burning 
ceases and the stars progress quickly as PG1159 onto the 
WD cooling track, while luminosity and mass-loss decrease 
and the wind reaches very high terminal velocities. 

A solid determination of the photospheric and wind pa- 
rameters, as well as of the chemical composition of CSPNe, 
is crucial to address questions concerning the stellar evolu- 
tion, the connection among the classes of H-deficient central 
stars, the wind driving mechanism and its effect on the ISM, 
and to understand the surrounding nebulae. Modelling the 
observed spectrum with stellar atmosphere codes is the most 
direct method for the derivation of stellar parameters, and 
constructing model grids covering the relevant range of pa- 
rameters is the best way to accomplish a systematic analysis, 
and constrain the major parameters. 

The majority of CSPNe spectral analyses in the past 
have been carried out in an ad-hoc fashion, that is, sets 
of models being calculated for each observed CSPN. These 
models vary in sophistication, with more or less ions be- 
ing included in the calculations, allowing or not for wind 



clumping, etc (see, for example, Leuenhagen et al. 1996 



Koesterke and Hamann1|1997b| |Leuenhagen and Hamann 
1998| |Koesterke et al.|1998||De Marco et al.|2001||Herald and 



Bianchi|2004a||Herald et aL|2005||Marcolino et al.|2007| ). To 

date, only a few Galactic CSPNe are well studied. With the 



sands of Milky Way post-AGB objects, including CSPNe 
and WDs, are being measured photometrically in the UV 
(Bianchi et al.||2005| [2007bl |2011a|b| ). Therefore, a general 



tool to facilitate determination of the stellar parameters will 
allow researchers to exploit the many-fold increase of sam- 
ples of these objects and to finally clarify the final stages 
of stellar evolution. We constructed a large grid of models 
covering the CSPN parameter space, which can be used to 
quickly constrain the main parameters and limit the need of 
'ad-hoc' computed models in each case. 

Several grids of stellar atmosphere models have been 
calculated to help analyse data for massive early type stars. 
Examples include: for main-sequence, giants, and supergiant 
stars, the "Kurucz" (in LTE and hydrostatic equilibrium) 
models |Kurucz|p9lj); for WN WR stars, the non-LTE 
grids of Hamann and Grafener (2004), which account for 
winds; for O-stars, the non-LTE, hydrostatic equilibrium 
models TLUSTY OSTAR 2002 ( |Lanz and Hubeny||2003l ), 
the non-LTE, accounting for expanding atmospheres, WM- 
basic grids of |Bianchi and Garcia (in prep. I and, in a more 



limited fashion, Pauldrach et al. ( 2001 ) ; for B-stars, the non- 
LTE, hydrostatic equilibrium TLUSTY BSTAR 2006 mod- 



els (|Lanz and Hubeny||2007| ). For CSPNe, Cle gg and Mid- 
|dlema ss (1987) computed a grid of non-LTE, plane-parallel 
model atmospheres, in hydrostatic equilibrium, with H and 



He opacity sources. Also, Rauch (2003) presented a grid of 
non-LTE, plane-parallel, line blanketed models calculated 
in hydrostatic equilibrium, useful for the analysis of H-rich 
CSPNe and their PNe. Only relatively recent advances in 
both computing power and stellar atmosphere codes made it 
possible to generate models that adequately reproduce most 
of the wind features in the spectra of CSPNe with winds. 
As a consequence, grids of non-LTE models, accounting for 
line blanketed expanding atmospheres and clumping, span- 
ning this important phase of post-AGB evolution do not yet 
exist. The grid presented here is the first step towards filling 
this gap. 

In this paper, we present a grid of synthetic spectra for 
the analysis of spectroscopic data from [WC] CSPNe hotter 
than 50000 K. The paper is organized as follows. Section [2] 
describes the stellar atmosphere code. We describe the grid 
in section [3] and present illustrative results in section [4] In 
section[5j we apply the grid to derive physical parameters of 
NGC 6905 in order to illustrate the use of the model grid. 
The paper finishes with conclusions in section [6] 



2 THE CODE 

[WC]-type CSPNe are very hot, with temperatures rang- 
ing from ~ 20000 to 200000 K and extended atmospheres, 
which expand reaching supersonic velocities. Line blanket- 
ing affects their atmospheric structure and emergent spec- 
tra. Not all CSPNe have dense winds, since the presence of 
intense mass loss, if driven by radiation pressure, is related 



to the closeness of the star to the Eddington limit (Paul- 
|drach et al.|1988"| . For example, as the PG1159 central stars 
approach the WD cooling track, their winds fade and their 
spectral lines are less and less conspicuous (see examples in 
|Werner et "aL"1|2010[ ). The winds of CSPNe are believed to 
be radiatively driven and as such, subjected to instabilities, 
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which are theoretically predicted to lead to the formation 



of clumping and shocks emitting soft X-rays (Owocki et al 



19881. Such effects were quantified observationally in the 



winds of massive O-type stars (|Bianchi and Garcia 2002 
Garcia and Bianchi||20fJ4} |Bianchi et al.||2009| . The grid of 



models was calculated using the non-LTE radiative transfer 
code CMFGEN ( jHillier and [ Miller|1998[|Herald et al.|2005| ), 
which solves iteratively for the level populations and radi- 
ation field, assuming radiative equilibrium, in a spherically 
symmetric expanding outflow. It accounts for wind clump- 
ing, soft X-rays, and line blanketing through a "super-level" 
approach. CMFGEN was originally developed to model the 
winds of massive WR stars and has been adapted for objects 
with weaker winds such as O-stars and CSPNe (as described 
by |Hillier et al.|2003| [Herald et al.|2005| ). Details about the 
code are given in the references above. Here, we briefly de- 
scribe the more important aspects. 

The radiative luminosity is taken to be constant 
throughout the wind and is given by 



(1) 



where R, and T* are defined at a Rosseland mean optical 
depth of 20. 

CMFGEN does not solve the dynamical equations of 
the wind, requiring the mass-loss rate and the velocity law 
to be supplied. Here, we have adopted a standard velocity 
law, where 



v(r) 



(2) 



with P — 1 which is the value usually adopted in modelling 
these objects, ro is roughly equal to R„. 

Currently, CMFGEN requires the density structure to 
be given. The procedure adopted here was to attach the 
wind velocity law to the underlying hydrostatic structure of 
the star (as described in Hillier and Lanz||2001 1, which we 
obtain using the TLUSTY code (Hubeny and Lanz 19951. 



TLUSTY computes plane-parallel, non-LTE atmospheres in 
hydrostatic equilibrium and requires surface gravity, tem- 
perature, and abundances as inputs. 

Instabilities in the stellar winds of hot stars are thought 
to generate density inhomogeneities, stochastically dis- 
tributed throughout the wind, which are usually referred 
to as clumping (see, for example, Moffat 2008). In CMF- 
GEN, the radiative transfer in an inhomogeneous medium 
is implemented assuming that the clumps are small com- 
pared to the mean free path of the photons and that the 
interclumping medium is void (this treatment of clumping 
is sometimes referred to as microclumping approximation). 
The degree of clumping in the wind is parametrized by the 
clumping filling factor /, with the density inside the clumps 
being a factor f^ 1 higher than the density of the homoge- 
neous wind model with same mass-loss rate. Since, accord- 
ing to this approximation, there is no material in between 
the clumps, the statistical equations are solved for the in- 
traclump medium. From the ion-level populations obtained, 
the non-LTE opacities and emissivities for the clump ma- 
terial are calculated. The transfer equation is then solved 
substituting the homogeneous-wind opacities and emissiv- 
ities with the ones calculated for the intraclump material 
multiplied by the filling factor. 

The main effect of treating clumps in the above ap- 



proximation is that it reduces the empirical mass-loss rate 
by a factor \J f~ x when using diagnostics that depend on the 
square of the density, such as emission lines and radio ther- 
mal emission. There is evidence of clumping in the winds 
of [WC] stars, but the quantification of the clumping fac- 
tor and its radial variation is made difficult by the weaker 
electron scattering wings of intense emission lines of these 
stars in comparison with their massive WR stars counter- 
parts. Despite that, Todt et al. ( 2008 1, based on the electron 



scattering wings of emission lines of three early-type [WC] 
stars, NGC 6751, NGC 5189, and NGC 1501, found evi- 
dence for clumping with /> 0.25. Further evidence for the 
presence of inhomogeneities in these winds comes from the 
work of |Grosdidier et aI7| ( 2000 , 200l|, who interpreted mov- 
ing features seen on the top of emission lines of the [WC] 
NGC 40 and BD +30 3639 as larger scale outflowing blobs. 
CMFGEN accounts for clumping using an exponential law, 



/(«) =/oo + (! -/oo) eX P(- U ( r )/«clump), 



(3) 



in which the radiative instabilities are supposed to be 
damped at low wind velocities and, as a consequence, the 
wind is assumed smooth until a certain wind velocity, v c i um p, 
is achieved. Beyond this velocity, the degree of clumping 
increases outwards until the clumping filling factor finally 
reaches its terminal value, f x , which, in this work, is as- 
sumed to be 0.1. In equation[3] f(v) and v(r) are the clump- 
ing filling factor and the velocity of the wind at a given 
radius, respectively. Here, u c i ump was taken to be 200 km 

CMFGEN accounts for line blanketing, which alters the 
ionization and temperature structure of the model atmo- 
sphere. In particular, the opacity by, e.g., iron lines redis- 
tributes the UV radiation field. Neglecting it would lead to 
an overestimation of the ionization of the atmosphere. The 
back warming and surface cooling effects of the line blanket- 
ing allow lower ions in the outer layers to coexist with the 
higher ones from the deeper layers. Inclusion, in the mod- 
els, of this phenomenon was noticed to better reproduce the 
spectra of massive WR stars ( Hamann and Grafener|2004 1 



In CMFGEN, the inclusion of the line blanketing is done 
through the superlevel approximation, in which several en- 
ergy levels are grouped into a smaller number of superlevels 
and all the levels within a superlevel share the same depar- 
ture coefficient. 

CMFGEN can also include, in the ionization calcula- 
tions, the effect of soft X-rays, which would be created in 
shocks distributed throughout the wind. These shocks are 
believed to originate from the instability of the line driv- 
ing mechanism. X-rays may considerably alter the ioniza- 
tion structure in the atmosphere, since they produce higher 
ionization stages of certain ions as, for example, O VI and 



N V (see, for example, 


Bianchi and Garcia|2002 


Garcia and| 


Bianchi 2004; Martins et al.)2005; Bianchi et al. 


|2007a| and 



references therein). In CMFGEN, the X-rays are taken into 
consideration in the calculations, by specifying the tempera- 
ture and velocity of the shocks. Also, a volume filling factor 
for the distribution of X-rays sources has to be given in or- 
der to set the level of the emission. The X-ray emissivity is 



taken from the X-ray code of Raymond and Smith (1977), 



which at present is available only for solar abundances. In 
section |4~H we describe the effect on the model spectra, of 
including soft X-rays in the calculations. 
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In order to calculate the detailed emerging spectra we 
adopted a microturbulence velocity which varies with depth 
such that 



-Rt — R* 



+ (v n 



v(r) 



(4) 



where « m in and f ma x are the minimum and maximum mi- 
croturbulence velocities. We assumed Vmin = 10 km s _1 and 
Wmax = 50 km s _1 for all the models. The models resolution 
varies across the wavelength range, with lines being bet- 
ter sampled than continuum regions. As an example, in the 
model B150.M70.v2500 (see Table [2), the central regions of 
the emission components of the O VI AA 1031.9, 1037.6 A 
and C IV AA 1548.2, 1550.8 A lines are sampled by 81 and 
46 flux points in an 1 A interval, respectively. 



3 THE GRID 

The grid presented here is intended to comprehensively cover 
the parameter space of [WC] CSPNe. The combination of 
parameters chosen for each of the grid models approximately 
follow the evolutionary calculations of Miller Bcrtolami and 



Althaus (20061 for CSPNe with final masses of 0.5, 0.6 and 



0.9 M© (which correspond to initial masses of 1.0, 3.1, and 
5.5 Mq, respectively), as shown in Fig. [I] The stellar param- 
eters L, R t , T t , loggr, and M» of the grid models are inside 
the range of values predicted for these stars. A chosen model 
temperature corresponds to a different model radius for each 
of the evolutionary tracks used, since the luminosity varies 
with mass on the constant luminosity track in the HRD. 
Also, for the models within one track, every temperature 
corresponds to a different radius in order to keep an almost 
constant luminosity. 

We have selected eight stellar temperature values rang- 
ing from 50000 to 200000 K and nine values of log g between 
4.0 and 7.0. For each combination of log g and T* seen in Fig. 
[I] models with different mass-loss rates and wind's termi- 
nal velocities were computed. The wind's terminal velocities 
and the mass-loss rates were chosen to cover the typical val- 
ues of these parameters found in the literature, taking into 
account that our models assume a clumped wind whereas 
some mass-loss rates were derived in the literature assum- 
ing a smooth wind. Additional models differing only in neon 



abundance were also computed (discussed in section 4.3 1 



The grid adds up to a total of 199 models. The list of mod- 
els is given in Tables [l] [2] and [3] The grid is available on- 
line at http://dolomiti.pha.jhu.edu/planetarynebulae.html 
There, the user will have access to the synthetic spectra and 
related documentation, along with plots comparing the dif- 
ferent models such as shown in Fig. [2] 

By following three tracks for different stellar mass val- 
ues, we covered different combinations of temperature and 
surface gravity /stellar radius. The [WC] CSPNe with tem- 
peratures above 50000 K, are characterized by the lack of 
strong absorption lines not affected by wind emission, which 
prevents the determination of their surface gravity through 
spectral analysis. When the distance is known, the derived 
stellar parameters can place the object on a specific evolu- 
tionary track, from which log g can be inferred. Otherwise, 
the analysis can only constrain the transformed radius value 
and mass and log g are not uniquely constrained. The trans- 
formed radius is defined as 



Uoo/2500 km s _1 



-i 2/3 



Af/10" 4 M yr- 1 



(5) 



and is a measure of how dense the wind is (smaller values 
of Rt translate into denser winds). Models with the same 
temperature, differing only in mass-loss rate and in stellar 
radius, but with the same transformed radius and the same 
wind's terminal velocity, are known to produce very similar 



wind features in the UV range ( Schmutz et al. 1989 Hamann 
et al. 1 1993 1 



Models having the same temperature, mass-loss rate 
and wind's terminal velocity, but from different mass tracks, 
differ in luminosity, radius, logg, and transformed radius. 
Different log g result in differences on the photospheric lines. 
Differing Rt affects the wind lines. Fig. [3] shows models 
A100.M67.V2000, B100.M67.V2000, and C100.M67.V2000 
that, despite having the same mass loss and wind's terminal 
velocity, have very different wind features as a result of their 
different transformed radii. These three models also differ in 
logg, but the lack of spectral features free from wind effects 
makes it impossible to distinguish among different log g val- 
ues. Surface gravity, however, influences the emission line 
profiles by affecting the underlying photospheric ones. 

It is believed that [WCL] and [WCE] stars form an evo- 
lutionary sequence due to their locations in the log T c ff-log g 
diagram. Carbon and helium are the main constituents of 
their atmospheres. Spectral analyses have, however, sug- 
gested different C:He ratios between the two subtypes 
(|Koesterke and Hamann||1997a|b|>. |De Marco et al.| ( |2001| ), 
Crowther et al.| ( |2003| >, and |Marcolino et al.| ( |2007| ), on the 
other hand, find no systematic discrepancies between the 
C:He ratios between [WCL] and [WCE] stars. Thus, since it 
is not yet clear if the two groups of [WC] CSPNe have dif- 
ferent carbon and helium abundance patterns, we adopted 
a constant C:He ratio for all the models in the grid. Typical 
measured values for elemental abundances are (by mass): 
X He = 0.33-0.80, X c =0.15-0.50 and X o =0.06-0.17 ( [Werner 
2001 [ |. Here, we adopted a C:He ratio of 0.45 : 0.43 and 
an oxygen abundance of Xo=0.08 by mass. Nitrogen abun- 
dances in these objects typically range from undetectable 
to ~ 2 per cent by mass, and we have adopted Xn=0.01 



(Koesterke and Hamann 1997b Leuenhagen and Hamann 



1998 1. We assumed a neon abundance of 2 per cent by mass, 



i.e., higher than the solar value by a factor of 11.5, since 
strong overabundances of this order have been reported for 



[WC] stars (Leuenhagen and Hamann 1998 Herald et al. 
2005[). Iron is expected to be depleted through s-process 
(jHerwig et al.| [2003j), as was indeed observed by |Stas ihska 
et al. (20041. Thus, the grid models have an iron abundance 



a factor of 100 lower than the solar value. Solar abundances 
were adopted for all other elements present in the models. 

The model calculations include many ionic species that 
have been previously neglected. The elements and ionic 
species considered in each model of the grid, along with the 
number of levels and superlevels used, can be found in the 
on-line documentation. The ionic species included in each 
model vary, since in many cases their number was limited 
to keep the models within a workable size, or due to un- 
availability of precise atomic data. The hotter models were 
calculated first and, as the temperature was decreased for 
the calculation of the cooler models, ionic species were in- 
troduced as needed based on the analysis of the ionization 
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fractions. Table[4]shows the species considered in the calcu- 
lation of the grid models. 

As a general rule, we considered it not necessary to in- 
clude a lower ionization stage when the ionization fraction 
of the lowest included ionization stage was < 1CT 2 . Tests 
were performed on representative models in order to assess 
the necessity of including lower ionization stages. As an ex- 
ample of the procedure, we give here more details about the 
C III, C IV, and C V ions. All models with temperatures 
up to 80000 K include C III. For temperatures of the order 
of 125000 K, C IV is much less abundant than C V; the 
bottom panel of Fig. [8] illustrates the typical behaviour. For 
some T„ =125000 K models showing relevant C IV ionization 
fractions, C III was included, but proved to be negligible and 
to have no important effect on the ionization fractions of C 
IV and C V. In models with T, = 100000 K, C IV may be- 
come comparable to, or more important than C V, at some 
radial distance from the star. In the T» =100000 K models 
where this distance is not very large, we added C III in the 
calculations. As shown in Fig. |1J the ionization fraction of 
C III is only ~ 10 -4 in the outer parts of the wind and even 
lower in the inner wind layers. Its inclusion caused no notice- 
able changes in the ionization fraction of the C IV ion, and 
was therefore considered not essential in this temperature 
regime. 



4 SPECTRAL DIAGNOSTICS OF M AND T, 

The preferred method of determining stellar temperatures 
is through the ionization balance of He and CNO elements. 
The observation of spectral lines from different ionization 
stages of the same element avoids the need for assuming 
abundance ratios. Consistence is achieved by the use of sev- 
eral diagnostic elements. However, in the spectra of [WC] 
CSPNe, not many elements show lines from different ioniza- 
tion stages. Therefore, the existence or absence of features 
in the spectra due to ions of different ionization potentials, 
along with the general appearance of the spectral lines, is 
also used as temperature diagnostic. 

In these stars, both photospheric and wind parameters 
affect the spectral features. Since several parameters can 
affect a spectral line, the identification of lines which are 
particularly sensitive to one parameter in a given regime 
is important for constraining its value. By examining the 
far-UV, UV and optical synthetic spectra produced for the 
grid, we identified spectral features sensitive to mass-loss 
rate and stellar temperature and selected those most useful 
for placing constraints on their values. These are discussed 
in sections 14. II and 14.21 



4.1 [WCE] types 

The far-UV, UV and optical spectra of [WCE] CSPNe are 
particularly poor in lines of multiple ionization stages of the 
same element. Only oxygen, neon and nitrogen show lines 
from multiple stages, yet only oxygen presents strong lines 
throughout the whole temperature range considered. 

Several other lines can be used as diagnostics of photo- 
spheric and wind parameters. The He II lines A 1640.4 and 
A 4685.7 A (the last one appears blended with a strong C IV 



line at A 4658.9 A, since in [WCE] stars, the wind's termi- 
nal velocities are high) are mostly sensitive to mass loss, and 
show little sensitivity to temperature in the regime between 
T* = 100000 and 165000 K. This happens because, for these 
high temperatures, helium is almost totally ionized. Exam- 
ples of line profiles are shown in Fig. [5j along with the helium 
ionization fractions, in which we see that the fraction of He 
II remains virtually unaltered in the hotter models. 

Among all the spectral lines in the studied range, Ne 
VII A 973.3 A shows the least sensitivity to mass loss, which 
makes it a good indicator of temperature within the range of 
physical parameters studied here. Its behaviour is shown on 
the top-left panel of Fig. [6] H erald et ah] ( |2005[ ) and |Bianchi| 



and Herald (2007) found, however, that it can be sensitive 



to the neon abundance in PC 1159 spectra, when it is not 
saturated. In section |4.3| we discuss the effect of the neon 
abundance in the grid models. 

For the O VI AA 1031.9, 1037.6 A doublet shown in 
Fig. [T] (top-left panel), both temperature and mass loss al- 
ter significantly the line profiles. The O VI A A 3811.4 and 
3834.2 A lines, seen in Fig. [7] (bottom-right panel), shows 
little variability with mass loss, within the range of values 
analysed, in the two cooler models. In the hotter models, 
on the other hand, it shows a behaviour opposite to that of 
all other strong lines in the range studied here, the line in- 
tensity decreasing with increasing mass loss. Also, the large 
variation of the line profile between models with tempera- 
tures of 125000 and 150000 K can help establish lower or 
higher limits for the temperature. In the bottom panels of 
Fig. [7] we show oxygen ionization fraction plots for differ- 
ent temperatures and mass-loss rates. We see that the O 
VI ion gains importance as we progress towards higher tem- 
peratures, which is reflected in the increasing intensity of 
the O VI lines in the hotter models. The mass-loss rate, on 
the other hand, has a more complex effect on these lines: 
it alters not only the density of the ion in the expanding 
atmosphere, but also changes its temperature structure, al- 
tering the ionization fractions. Denser winds will be cooler, 
with lower ionic species growing in importance. For mas- 
sive stars, the O VI AA 1031.9, 1037.6 A lines were found 
to strongly depend on X-rays and to have a hard threshold 
with mass loss, while other lines strongly depend on clump- 
ing (Bianchi and Garcia 2002). The effect of X-rays on the 



far-UV O VI doublet will be discussed in section EL4l 

Other temperature indicators are the neon and oxygen 
features shown in Figs. [6] (top-right panel) and[7] (central-left 
panel), respectively. These features are due to O IV (with 
lines at AA 3403.6, 3411.7, and 3413.6 A), O VI (at A 3433.3 
A) and Ne VI (with contamination by Ne V lines). The Ne 
VI feature can only be seen in hotter models (T* 150000 
K). The O IV feature appears only in the two cooler ones (for 
T* = 100000 and 125000 K), while it is substituted by the 
O VI line as the temperature increases. For the right combi- 
nation of mass loss and temperature, these oxygen lines can 
also coexist. Thus, these features can be used to constrain 
the temperature regime of the object being studied. Also 
seen in the central-left panel of Fig. [7|are the Ne V A 3313.7 
A and N IV AA 3478.7, 3483.0, and 3484.9 A lines, which 
can help constrain neon and nitrogen abundances. 

The strong C IV doublet at AA 1548.2, 1550.8 A shows 
sensitivity to both mass loss and temperature, as can be seen 
in the central-left panel of Fig. [8] Nevertheless, once either 
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mass loss or temperature is established through other indi- 
cators, it can help constraining the other parameter. Also, 
this line is an important diagnostic of wind's terminal ve- 
locity, since this ion is abundant also in the outer parts of 
the wind down to temperatures of 50000 K. 



tracting that of the absorption component from that of the 
emission. These two lines are very sensitive to mass loss in 
the [WCE] and [WCL] parameter regimes, respectively, and 
thus can be used to identify the interval of transformed ra- 
dius and temperature that matches the data. 



4.2 [WCL] types 

The winds of [WCL] CSPNe are characterized by much lower 
terminal velocities than the ones of [WCE] stars, which re- 
sults in narrower spectral lines. Also, these cooler CSPNe 
seem to show a wider range in mass-loss, by comparing re- 
sults of spectral analyses found in the literature. We took 
these facts into consideration when calculating the grid. 

These objects present very complex spectra, and a 
higher number of spectral lines from elements with multiple 
ionization stages, thus making the stellar temperature easier 
to constrain. Several strong C III and C IV lines appear in 
the grid models for these stars, along with multiple ioniza- 
tion stages of neon, nitrogen, oxygen, phosphorus, sulphur, 
and helium. 

We selected some spectral features useful for distin- 
guishing among temperatures and mass-loss rate values, 
within the regime considered. The top-left panel of Fig. [8] 
shows a section of the optical spectra where strong lines 
from both C III (A 5695.9 A) and C IV (AA 5801.3, 5812.0, 
and 5865.9 A) can be seen, along with a He I line (A 5875.7 
A), making this spectral region an important diagnostic of 
stellar temperature. The O IV AA 3403.6, 3411.7, and 3413.6 
A spectral features shown in the central-right panel of Fig. 
[7] is likewise an excellent indicator of stellar temperature 
between T» = 50000 and 80000 K, showing very little sensi- 
tivity to mass-loss rate in the interval studied. Also shown 
is the N IV A 3478.7 A line. 

In the late-[WC] type, the O VI doublet AA 1031.9, 
1037.6 A, shown in Fig. 171 top-right panel, is still present, 
albeit much weaker and showing contamination by other 
spectral lines. Another useful diagnostic line for the [WCL] 
subclass is the C III A 1908.7 A line shown in the top-right 
panel of Fig. [8] It is much more sensitive to mass loss than 
to temperature in the range shown, thus making it an im- 
portant mass loss discriminator. 

Other lines such as He II A 1640.4 A (shown on the 
top-right panel of Fig. He II A 4685.7 A and C IV A 
4658.9 A (central-right panel of Fig. [HJ, O III AA 1150.9 
and 1153.8 A (bottom-left panel of Fig. [7J also showing a 
Ne V line at A 1146.1 A), C IV ~ A 1169 and C III ~ 
A 1176 A (both of them are blends of several lines) and 
C III A 2296.9 A (Fig. [£] central-right and bottom panels) 
show high sensitivity to both mass loss and temperature. 
Nevertheless, once temperature or mass loss are constrained 
by other indicators, they can be an important diagnostic of 
the remaining parameter. 

In order to further visualise the dependence of spectral 
diagnostics on the photospheric and wind parameters, we 
made contour plots of the equivalent widths of two spec- 
tral lines discussed in this section, having stellar tempera- 
ture and transformed radius on the axes. Contour plots of 
the He II A 1640.4 A line, and the C III A 1908.7 A line 
are shown in Fig. |9j where the contours indicate the equiv- 
alent widths of the lines. In the case of P-Cygni profiles, 
CMFGEN calculates the equivalent widths of a line by sub- 



4.3 Effect of Ne abundance 

A neon abundance of the order of 2 per cent by mass is ex- 
pected from evolutionary models for the [WC] stars, which 
present in their surfaces the abundance pattern of the re- 
gion between the hydrogen and helium burning shells of the 
precursor AGB star. In AGB stars, neon is produced from ni- 
trogen through the 14 N(a, 7) 18 F(/3 + ) ls O(a, 7) 22 Ne chain in 
the helium burning shell and then dredged up to the convec- 
tive intershell (see, for example, Werner and Herwig 2006). 
Thus, the determination of neon abundances is important 
to constrain evolutionary models, in particular, the dredge- 
up episodes. Also, the presence of several ionization stages 
of neon in the winds of [WCE] stars is potentially a use- 
ful tool for the determination of stellar temperature. In Fig. 
[6] we show several neon features present in the models in 
different temperature and mass loss regimes, for different 
values of neon abundance, ranging from solar to 11.5 times 
super-solar (which is adopted for the models in the grid). 
On the top-left panel, the Ne VII A 973.3 A line seen in the 
grid models for [WCE] stars shows little sensitivity to neon 
abundance in the models with T» = 100000 and T* = 125000 
K. For higher temperatures, this line starts to show a higher 
sensitivity to Ne abundances and can potentially help con- 
strain it. Also, if used together with neon lines from different 
ionization stages, it can contribute to the determination of 
the stellar temperature, although it also depends on mass- 
loss rate. 

In the central-right and left panels of Fig. [6] we show 
the Ne III lines, which appear in the synthetic spectra 
for the [WCL] sub-class at A 2553.4 A and at AA 2677.9, 
2678.6 A, respectively. Both lines present high sensitivity to 
neon abundance, temperature, and mass-loss rate, and are, 
therefore, useful to constrain the neon abundance when the 
photospheric and wind parameters have been determined 
through other indicators. They also help the determination 
of stellar temperature when used in conjunction with other 
neon lines of different ionization stages. 



4.4 Effect of X-rays 

X-rays emission of the order Lx ~ 10 -7 L„ observed in mas- 
sive O and B-type stars ( Chlebowski et al.|1989 Chlebowski 



and Garmany|1 991 Evan s et al.|2003 1 is commonly believed 
to originate from shocks in the stellar wind, which form due 
to the unstable nature of the radiative force on spectral lines 



( |Lucy and White||1980[|Owocki et aL| [l988). Such a flux of 
X-ray photons has been invoked in order to reproduce the 
strong O VI AA 1031.9, 1037.6 A doublet in early O-type 
stars and to achieve consistency among observed N V, Si 



IV and C IV diagnostics ( Garcia and Bianchi|2004 Bianchi 
et al.|2009 i. 



Since the fast winds of CSPNe are also believed to be 
subjected to the radiation pressure in spectral lines, they 
are expected to be clumped and may exhibit X-ray fluxes. 
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Guerrero et al. (20011 measured an X-ray luminosity of 



10~ L* for the central star of NGC 6543, consistent with 
that predicted from shocks in the stellar wind. However, the 
emission could also be explained by coronal activity of an 
undetected companion. 

We examined the effect of soft X-ray fluxes in the syn- 
thetic spectra computed for our grid. We re-computed mod- 
els with T» = 65000, 100000, and 165000 K including soft 
X-ray fluxes into the calculations. We assumed shock tem- 
peratures of 5 x 10 6 and 3 x 10 6 K, starting at velocities 
of 300 and 400 km s -1 , respectively. In the model with 
T* = 165000 K, only an extremely high X-ray flux yielding 
an observed X-ray luminosity of ~ 10~ 2 L* (for a mass-loss 
rate of 1.0 x 10 -7 Mq yr -1 ), caused a relevant effect on the 
spectra, mainly in the O VI doublet AA 1031.9, 1037.6 A, as 
shown in Fig. |10| In models with progressively lower temper- 
atures, we found that lower X-ray luminosities modify the 
O VI doublet profile. For the T* = 100000 K model, a rea- 
sonable X-ray luminosity of ~ 10~ 7 L* (for M = 1.0 X 10 -7 , 
M = 2.0 x 10" 7 , and M = 3.0 x 10" 7 M yr" 1 ) significantly 
alters the O VI doublet (but not other spectral features), as 
shown in Fig. [10] At a temperature of T„ = 65000 K the 
convergence of models with X-ray fluxes proved to be very 
hard to reach. Nevertheless, we were able to obtain a con- 
verged model for a mass-loss rate of M = 5.0 x 10 -8 Mq 
yr -1 with an X-ray luminosity of 4 x 10~ 10 L„. Again, the 
only spectral feature affected was the far-UV O VI doublet, 
which is also shown in Fig. |10| 



5 SPECTRAL ANALYSIS OF THE CENTRAL 
STAR OF NGC 6905 

To illustrate the usefulness of the model grid, we apply it 
here to derive the physical parameters of the [WCE] central 
star of NGC 6905. We compare our grid models with data 
(Table [5| available at the Multimission Archive at STScI 
(MAST) from FUSE and HST spectrographs. FUSE cov- 
ers the wavelength range 905-1187 A with a resolving power 
of < 20000. We also used HST-STIS spectra obtained with 
G140L (R ~ 1200) and G230L (R ~ 750) gratings with wave- 
length coverage of 1150-1736 A and 1570-3180 A, respec- 
tively. In all figures showing the observed spectra of NGC 
6905, the resolution of the models, in the STIS wavelength 
range, is matched to the resolution of the observations by 
convolving the synthetic spectra with the instrument line 
spread functions, while for the spectral region covered by 
FUSE, which is full of interstellar hydrogen absorptions, the 
models and observations were convolved with a Gaussian of 
FWHM = 0.1 A to facilitate the visualization. In Figs. 
|11|20| the interstellar absorption was omitted from the syn- 
thetic spectra to facilitate visualization and will be further 
discussed below. The line identifications in these figures only 
mark the strongest lines in each model. 

We obtained an initial estimate of the wind's terminal 
velocity by measuring the blue edge velocity, « e d ge , from 
the strong P-Cygni profile of the C IV A 1548.2 A line. We 
derived = 2170 km s~ x , correcting v e d ge for an assumed 
turbulence velocity, u tU rb, of 10 per cent v^, (Lamers and 
Cassinelli|1999 l. 



fedge 2 X Vturb- 



(6) 



This estimate is confirmed by comparing the observed spec- 
tra with grid models computed for wind's terminal velocities 
Woo = 2500, 2000 and 1500 km s _1 . An example is shown in 
Fig. EH 

We then used the grid of models with the appropriate 
wind's terminal velocity to constrain the temperature and 
the transformed radius values of the central star of NGC 
6905. Based on these results, we then perform a more de- 
tailed analysis, varying additional parameters not explored 
by the grid. 



5.1 Using the grid to constrain temperature and 
transformed radius 

The central star of NGC 6905 is a [WCE] , therefore its spec- 
tra are poor in lines of multiple ions of the same element. 
Among the features identified in the spectra, only oxygen 
shows strong lines of multiple ionization stages. Several O 
VI lines are present (AA 1031.9, 1037.6 A, AA 1080.8, 1081.6 
A, AA 1122.5, 1124.8, 1124.9 A, A 1291.9 A, A 1996.1 A, A 
2070.4 A, A 2082.04 A and A 2431.409 A), along with two 
strong O V lines at A 1371.3 A and AA 2781.0, 2787.0 A. 
However, as will be discussed below, these two O V lines are 
only matched by grid models incompatible with other spec- 
tral diagnostics and both show sensitivity to the inclusion 
of heavy elements in the calculations, while other lines do 
not. Therefore, they are more useful as indirect indicator of 
abundance and not to uniquely constrain temperature. Sev- 
eral C IV lines are also observed and used as temperature 
diagnostic. 

As a first step, we compare the observed spectra of the 
central star of NGC 6905 with our grid models with a wind's 
terminal velocity of 2000 km s _1 , the closest value available 
in the grid to the one inferred from the C IV A 1548.2 A 
P-Cygni absorption edge. The comparison rules out mod- 
els with T* = 125000 K, which show very weak O VI AA 
1031.9, 1037.6 A, of about half the observed intensity, for 
all the mass-loss rate values in our grid, besides a too weak 
Ne VII A 973.3 A line, which, as shown in Fig. [6] at this 
temperature, shows little sensitivity to mass loss and neon 
abundance. The carbon, helium and other O VI lines can 
be simultaneously matched by models of this temperature. 
The OVA 1371.3 A is too strong in T„ = 125000 K models 
and shows a deep absorption component which is not ob- 
served. As we will show below, this line is sensitive to the 
inclusion of heavier ions in the models and their abundances 
and tends to become stronger as new ions are added to the 
models. Models with a temperature of 200000 K can also be 
ruled out since, for the mass-loss rate values that allow rea- 
sonable fits of the carbon and helium lines, both O V lines 
are absent and very strong unobserved Ne VIII features ap- 
pear in the synthetic spectra, which also show weak O VI 
lines. The models with temperatures of T* = 150000 and 
165000 K better agree with the observed spectra and are 
shown in Figs. [12] and |13| respectively, for different values 
of transformed radius. Among the T» = 150000 K models, 
the B150.M70.V2000 one with M = 10" 7 M yr" 1 , and 
Rt = 10.72 Rq fits best all the observed diagnostics, except 
for the O V AA 2781.0, 2787.0 A and A 1371.3 A lines, which 
are only fitted by models of this T» , by adopting a mass-loss 
rate twice and 1.3 times higher, respectively. Such a higher 
mass loss, however, worsens the fit of all the C IV, He II and 
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O VI lines, except for the O VI AA 1031.9, 1037.6 A doublet, 
where the intensity of the second component is improved, as 
shown in Fig. |12| The relative strength of this line, however, 
is also very sensitive to clumping, since emission lines are af- 
fected by it, but not resonance ones ( |Oskinova et al.|2007[ ), 
as well as to soft X-rays, and other factors. Among models 
with temperature of T* = 165000 K, the C165.M65.V2000 
one, with a mass-loss rate of M = 10 -7 Mq yr -1 , and a 
transformed radius of _R t = 9.03 fits well the observed 
spectra, again except for the O V lines which are even more 
discrepant than in the T„ = 150000 K models, as can be 
seen in Fig. |13| At any of these two temperatures, the fits 
of the O V lines require lower values of transformed radius, 
which are incompatible with the other diagnostic lines. This 
effect is worse in the T* = 165000 K models than it is in the 
T» = 150000 K ones, from what we conclude that, among 
the grid models, the B150.M70.V2000 one gives the best fit. 

Below, we extend the analysis beyond the model grid 
comparison and address in depth not only the O V lines dis- 
crepancy but also other disagreements between the observa- 
tions and the best-fitting model, and refine the parameters. 
We compare the observations with synthetic spectra calcu- 
lated assuming different values of turbulence velocity, show 
the effect of different neon, oxygen and argon abundances, 
put constraints on the nitrogen and iron abundances, and 
study the impact of adding heavier ions in the synthetic 
spectra, which we find to affect the ionization fractions of 
oxygen and to improve the O V lines fit. 



5.2 Extending the analysis 



In Fig. 
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the O VI AA 1031.9, 1037.6 A doublet and the ab- 
sorption profile of C IV AA 1548.2, 1550.8 A are narrower m 
the Woo = 2000 km s _1 models than in the observations, sug- 
gesting that the turbulence velocity in NGC 6905 is higher 
than that adopted in the calculation of the model grid. Fig. 
|14| shows that higher values of turbulence velocity, in fact, 
improve the match of the O VI and C IV profiles. 

Another difference between the best-fitting grid model 
and observed spectra is the absence, in the observed spectra, 
of N V AA 1238.8, 1242.8 A seen in all the models. This is due 
to a lower nitrogen abundance in this star than we assumed 
for the grid. Fig. [15] shows models with different nitrogen 
abundances and T» = 150000 K, M = 10" 7 Mq yr" 1 , and 
Woo = 2000 km s _1 . Nitrogen abundances of An < 5.5 x 10 -4 
better match the spectral region shown in Fig. |15| while 
other spectral lines show no significant differences. 

Several neon features which appear in the synthetic 
spectra are not observed, indicating that this object may 
have a lower Ne abundance. Fig. [16] shows models with dif- 
ferent neon abundances and T, = 150000 K, M = 10~ 7 M Q 
yr -1 , and Voo = 2000 km s _1 . The neon abundance in the 
spectra of the central star of NGC 6905 is difficult to con- 
strain due to the lack of strong neon lines other than Ne VII 
A 973.3 A, and due to the mismatch between Ne VI and Ne 
V features present in the models in t he r egion between 1700- 
2300 . 



16 



we see that none of 



and the observations. In Fig. 
the neon abundances adopted improves the fit of the Ne VI 
complex AA 2213.1, 2229.1 A. Also, models with higher neon 
abundances produce some features which are not seen in the 
observed spectra such as a Ne VII line at A 2163.4 A and Ne 
VI lines at A A 2641.1, 2687.4 A. On the other hand, the lower 



neon abundances produce very weak Ne VII A 973.3 A and 
Ne VI AA 997.4, 999.6 A lines. The fact that some diagnostics 
are improved by a higher neon content while others are wors- 
ened seems to point to an incorrect temperature. We did not 
use the neon lines as a temperature diagnostic as they are, 
also in the T* = 165000 K models, incompatible with other 
diagnostics. However, we have observed (not shown here) 
that models with temperatures of T» = 200000 K better re- 
produce several of the neon structures that are incompatible 
in cooler models. We consider the present analysis inconclu- 
sive regarding the neon abundance. An interesting outcome 
of the comparison among models with different neon abun- 
dances is a slight enhancement of the OVA 1371.3 A and 
AA 2781.0, 2787.0 A lines as neon content in the models 
increases, which is the result of an increase in the relative 
fraction of the O V ion, as can be seen in the bottom panel 
of Fig. [16] 

In an attempt to reproduce the strong O V lines simul- 
taneously with the other diagnostics, we investigated the 
effect of higher oxygen abundances, which is shown in Fig. 
|17| For the grid models we assumed an oxygen mass fraction 
of Ao = 0.08. A higher oxygen abundance increases the in- 
tensity of the O V lines in the synthetic spectra, but worsens 
the fit of most of the O VI lines. Besides, even an oxygen 
mass fraction of Xo = 0.20 is not enough to reproduce the 
strong observed O V AA 2781.0, 2787.0 A line. The O VI 
AA 1031.9, 1037.6 A doublet does not show any significant 
change from the increase in the oxygen abundance. 

Because of the effect seen in the O V lines due to an 
increase in neon abundance, we also experimented with the 
inclusion of new ions in the models in an attempt to repro- 
duce the O V lines observed intensities. A very complete 
model including Ni VII, Ni VIII, Ni IX, Co VII, Co VIII, 
Co IX, Ca VI, Ca VII, Ca VIII, Ca IX, Ca X, Ar VI, Ar VII, 
Ar VIII, Si V, Si VI, Mg V, Mg VI, Mg VII, Na VI, Na VII, 
Na VIII, Na IX in addition to the ions already present in the 
previous model was calculated, assuming solar abundances 
for all the additional elements. As a result, the intensity of 
the O V A 1371.3 A and AA 2781.0, 2787.0 A lines in the 
synthetic spectra increased as shown in Fig. |18| while other 
major features which we used to constrain stellar parame- 
ters were unaffected. This is a result of the increase of the O 
V ion fraction. In particular, the magnesium ions followed 
by the sodium ones produced most of the effect, while the 
cobalt and nickel ions cause smaller, but noticeable changes 
of the O V line strengths. The rest of the spectra show minor 
differences, mainly due to calcium and argon lines. 

We further explore this model by altering the abun- 
dances considered for some elements and show the effect 
on the spectra of diffe rent abundances of the elements ar- 
gon and iron. Fig. 19 shows the Ar VII A 1063.55 A line. 



This line was first identified as a photospheric absorption 



by Werner et al. (20071 in CSPNe and WDs, for which they 



found roughly solar abundances, in line with nucleosynthe- 
sis calculations for AGB stars. We find a model with ten 
times the solar abundance to better reproduce the observed 
feature in this spectral region, confirming the results of |Her-| 
|ald and"Bi anchi ( 2009]) who identified this same line in the 



20 



FUSE spectra of some of the hottest [WC] CSPNe. Fig 
shows the effects of different iron abundances in the syn- 
thetic spectra. In a recent work, | Werner et al.| ( |2010[ ) iden- 
tified Fe X AA 979.3, 1022.9 A lines in the FUSE spectra 
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of the hot PG1159 stars RX J2117.1+3412, K 1?16, Long- 
more 4, NGC 246, H1504+65, finding a solar iron abundance 
for these stars. In NGC 6905, we find that iron abundances 
higher than 0.3 times the solar one produce some unobserved 
features, especially an absorption on the Ne VII A 973.3 A 
P-Cygni profile, and can be ruled out. 



5.2.1 Interstellar absorption 

Figs. [2l] and [22] show our best-fitting model, which in- 
cludes the additional heavy ions and an argon abundance 
of ten times the solar value, calculated adopting a turbu- 
lence velocity Wturb = 150 km s" 1 . We applied the inter- 
stellar absorption due to atomic and molecular hydrogen, 
treated as described in |Herald and Bianchi (2002, 20 04a|b[ ), 
to the synthetic spectra. We assumed that the tempera- 
ture of the interstellar gas is 100 K, a typical value for 
the ISM. From fitting the Ly Q A 1215 A H I absorption 
profile, we derived a neutral hydrogen column density of 
20.7 < log N (HI) < 21.1 (where N is given in units of 
cm -2 ). Using the numerous H2 absorption lines in the FUSE 
range, we constrained the molecular hydrogen column den- 
sity to 19.3 < log N(Hz) < 19.7 (where N is given in units 
of cm -2 ). 

After the stellar parameters have been constrained 
based on line diagnostics, the observed slope of the spec- 
tra was compared to the best-fitting model reddened us- 
ing different values of the colour excess, as shown in Fig. 
23] The reddening of|Cardelli et al.| (|1989} with Ry = 3.1 



also lower value for the distance, a radius of 7?* = 0.09 Rq 
is obtained. A model with this radius and T* = 150000 K 



was adopted. The STIS G140L and G230L spectral regions 
are well matched using a colour excess of -Eb-v = 0.17 
mag, while a lower reddening is required in order to fit 
the FUSE region. Thus, we conclude on a colour excess of 
Eb-v = 0.17 ± 0.05 mag. The continuum comparison sug- 
gests that the interstellar extinction may be less steep in 
the UV than the extinction curve adopted. Model contin- 
uum accuracy has not yet been reportedly tested versus ob- 
servations and may be addressed with future work on an 
expanded sample. 
According 



to 



Bohlin et al. (19781, the neutral hy- 



drogen column density and the colour excess relate as 
(N (HI) /Eb-v) = 4.8 x 10 21 atoms cm" 2 mag" 1 . An 
Eb-v = 0.17 mag implies log N (HI) = 20.9 (where N is 
given in units of cm" 2 ), which is within the range derived 
by us from the Ly a A 1215 A absorption. 



would be in a lower mass evolutionary track ( Miller Berto- 
|lami and~A lthaus 2006) of 0.57 M , as seen in Fig. 24 where 
the filled circle is the grid model calculated by us and the 
open circle is assuming the literature distance of 1.7 kpc. 
By keeping the transformed radius constant, the radius of 
R, = 0.09 R0 implies a mass-loss rate of M — 6.2 x 10" 8 
M Q yr" 1 . 

The central star of NGC 6905 has been previously stud- 
ied in the literature through spectral analyses in optical, 
ultraviolet and far-ultraviolet regimes. It has been recently 



analysed by Marcolino et al. (20071, using low resolution 



optical spectra, low resolution IUE spectra and FUSE spec- 
tra in the 1000 — 1175 A interval, which excluded the Ne 
VII A 973.3 A line from the analysis. They also made use of 
the code CMFGEN to perform the spectral analysis and 
found parameters close to the values determined by us: 
T» ~ 150000 K, M ~ 7 x 10" 8 M yr" 1 , R t = 10.5 Rq, 
and Vex, = 1890 km s" 1 , considering a clumping filling fac- 
tor of /= 0.1, the same value adopted by us. In our analy- 
sis, we included many ionic species not considered in theirs, 
which allowed us to achieve a good fit of the OVA 1371.3 
A line and improve the fit of the O V AA 2781.0, 2787.0 
A line. This example proves the usefulness of the grid for 
constraining the main stellar parameters, and our extended 
analysis resolved the so far inconsistent diagnostic of the O 
V lines. 



5.2.2 Results from the analysis 



The parameters of our best-fitting grid model are T, = 
150000 K, M = 1 X 10" 7 M Q yr" 1 , R t = 10.7 R , R, = 0.12 
Rq, and v x = 2000 km s~ . The value of the stellar ra- 
dius follows the evolutionary track of |Miller Bertolami and| 
Althaus| p006l for 0.6 M CSPNe and a temperature of 



150000 K. Nevertheless, as discussed previously, two mod- 
els sharing R t , i>oo , and T* , but with different values of _R» 
and M will have similar spectral features. As a consequence, 
different combinations of parameters result in equally good 
fits to the observations. In the literature, the distances de- 



rived range between 1.7 and 2.6 kpc ( 


Maciell 1984 


Cahn 


et al.||1992| |van de Steene and Zijlstra 


|1994| (Zhangl 


1995 


Stanghellini et al. 


2008 


). If we adopt the most recent and 
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6 SUMMARY AND CONCLUSIONS 

We have used the state-of-the-art stellar atmosphere code 
CMFGEN to build a comprehensive grid of synthetic spectra 
covering the range of stellar parameters appropriate for H- 
deficient CSPNe of temperatures above 50000 K. The mod- 
els account for line-blanketing, non-LTE expanding atmo- 
spheres (where the velocity law is fixed and the mass-loss 
rate is a free parameter), and wind clumping. Many ionic 
species neglected in previous models were included (Table 
|4| . We provide a uniform set of models that allows system- 
atic analysis of far-UV, UV, optical and IR observed spectra 
in order to constrain photospheric and wind parameters. We 
examined the far-UV, UV and optical synthetic spectra of 
the grid and selected spectral features best suited for con- 
straining values of mass-loss rate and stellar temperature. 

We also analysed the effect of soft X-rays on the syn- 
thetic spectra. We found that, at T» = 165000 K, an unlikely 
high X-ray luminosity was necessary in order to cause a rel- 
evant effect in the spectra. At T„ = 100000 K, a reasonable 
value of L x strongly affects the O VI AA 1031.9, 1037.6 A. 
For the T, = 65000 K model, an even lower value of X-ray 
luminosity is sufficient to affect the O VI doublet. 

As a sample application of the grid, we used it to analyse 
the UV spectra of the central star of NGC 6905 showing that 
the grid can be used to constrain the photospheric and wind 
parameters to a good extent and to greatly facilitate more 
detailed analysis. 

For NGC 6905, the temperature of our best-fitting 
model is also a lower limit, since any lower temperature 
model would produce a too strong OVA 1371.3 A line 
which, as we have shown, is also affected by the inclusion 
of new ions. A slightly higher temperature for this star im- 
proves the fits of the far-UV O VI doublet and Ne VII line, 
while the other diagnostics remain acceptable. Thus, we con- 
strain the temperature of the central star of NGC 6905 as 
being between 150000 and 165000 K, a 10 per cent overall 
uncertainty. 

We extended the study beyond the grid of models by 
considering parameters that were not explored in the grid. 
We found higher values of the turbulence velocity to im- 
prove the fit of O VI AA 1031.9, 1037.6 A and C IV AA 
1548.2, 1550.8 A lines. We showed the effects of different 
abundances of neon, oxygen and argon and we were able to 
put constraints on the nitrogen and iron abundances, both 
subsolar. In this object, a lower value of nitrogen abundance 
than the one chosen for the grid was necessary to fit the 
observed spectra. For all plausible temperatures, we were 
unable to adequately fit all neon diagnostics simultaneously 
with a single neon abundance. 

By adding heavier ions to the best-fitting grid model, 
we were able, for the first time, to reproduce the observed 
intensity of the OVA 1371.3 A line and to improve the fit 
of the O V AA 2781.0, 2787.0 A line. Finally, by fitting the 
several H2 absorption lines present in the FUSE range and 
the Ly Q absorption profile, we constrained the molecular 
and neutral hydrogen column densities and, by comparing 
the slope of the observed spectra to our best-fitting model, 
we derived the colour excess, which is in agreement with the 
value expected from the H I column density measured. 

The grid of synthetic spectra presented here is avail- 



able on-line, at http://dolomiti.pha.jhu.edu/ planetaryneb 



ulae.html. Additional models will be added to the grid in 
the future, covering the temperature regime below 50000 K, 
as well as new grids for different surface gravity regimes, 
suitable for the study of PG1159 stars that present wind 
signatures in their spectrum. 
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Figure 1. Grid points in the log T*-log g plane. The evolutionary tracks of Miller Bertolami and Althaus ( 2006 ) are also shown (continuous 
lines) for CSPNe with 0.5, 0.6, and 0.9 Mq, from right to left. Each point shown is identified by a label which corresponds to a group 
of models with the same temperature, surface gravity, and radius, but with a range of mass-loss rates and wind's terminal velocities, as 
can be seen in Tables [l] [2] and [3] 
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Track B: M = 0.6M„: V =2000 km s*~ 



Track B: M = 0.6M o ; v.= 2000 km s _I 




Figure 2. Examples of plots comparing different grid models, from the far-UV to the optical range. Here and in subsequent figures, 
mass-loss rates are given in Mq yr _1 . Similar plots for all grid models are available on-line. 



R, = 6.40R S ; logg = 6.0 (A100.M67.V2000) 
R,= lD.67R a ; logg = 5.3 (B100.M67.V2000) 
R,=a3.42R a ; logg=4.8 (C 100.M67.V2000) 
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Figure 3. Models A100.M67.V2000, B100.M67.V2000, and Cf 00.M67.V2000, with T = fOO kK, M = 2x 10" 7 Mq yr" 1 and Voo = 2000 
km s — 1 , differ in R, Rt, and logg. Wind features vary greatly due to the different wind densities. 
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Figure 4. The C ionization fractions of two T*=100000 K models, with identical stellar parameters, and the C III ion included (dashed 
lines) or not included (solid lines) in the calculations. 
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Figure 5. Model profiles of the He II lines A 1640.4 A (top panels) and A 4685.7 A (central panels) at different temperatures (also radii) 
and mass-loss rates. The bottom panel displays the ionization fractions of helium, where continuous lines indicate models with M = 10~ 7 
M yr" 1 , and dashed lines indicate M = 10~ 6 M Q yr" 1 if T* < 80 kK, or M = 3 X 10~ 7 M yr" 1 for models with T» > 100 kK. All 
the models shown belong to track B (M, = 0.6 Mq) and have Hoo = 2500 km s — 1 . 
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Figure 6. Ne lines from different ionization stages and their variation with temperature (also radius), mass-loss rate and neon abundance, 
fn the top-left panel, for clarity, the models with Xjmc = 0.01 were omitted when the effect of Ne abundance is small. The bottom panel 
shows the ionization fractions of neon (for models with a neon mass fraction of X-^ c = 0.02); continuous lines indicate models with 
M = 10~ 7 Mq yr _1 , and dashed lines indicate M = 10~ 6 M yr~ x for models with T, < 80 kK, or M = 3 X 1CT 7 Mq yr" 1 for models 
with T t > 100 kK. All the models shown are from track B (M„ = 0.6 M Q ) and have i) M = 2500 km s _1 . 
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Figure 7. Synthetic line profiles from different ionization stages of oxygen in various models. The bottom panel shows the ionization 
fractions of oxygen; continuous lines indicate models with M = 10 — 7 Mq yr — 1 and dashed lines indicate M = 10~ B Mq yr -1 if T t ^ 80 
kK, or M = 3 X 10" 7 Mq yr _1 for models with T* ^ 100 kK. All the models shown arc from track B (M, = 0.6 Mq) and have Voo = 2500 
km s~ - 1 . 
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Figure 8. Profiles of carbon lines of different ionization stages from models with different mass-loss rates and temperatures (also radii). 
The bottom panel shows the ionization fractions of carbon; continuous lines indicate models with M = 10 — 7 Mq yr _1 and dashed lines 
indicate M = 1CT 6 M yr" 1 if T» ^ 80 kK, or M = 3 X 10~ 7 M yr" 1 for models with T» > 100 kK. All the models shown are from 
track B (A/„ = 0.6 M ) and have v x = 2500 km s _1 . 



A New Grid of Synthetic Spectra for [WC] CSPNe 19 





Figure 10. Comparison between models with (dash-dotted line) and without (continuous line) X-rays. The O VI AA 1031.9, 1037.6 A 
doublet is shown due to its sensitivity to X-ray ionization (see text). For the T* = 165 kK model, an extreme X-ray luminosity of ~ 10~ 2 
L* is shown, for the T* = 100 kK models, we show L x ~ 10~ 7 L», and for the T* = 65 kK model, L x = 4 X 10~ 10 L» 
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Figure 11. The two lines showing the strongest observed (continuous black line) P-Cygni profiles in the spectra of NGC 6905, O VI 
AA 1031.9, 1037.6 A in the FUSE range and C IV AA 1548.2, 1550.8 A in the STIS G140L range, are shown, as well as models computed 
with three values of v^x, from our model grid. All other parameters in the models are identical and correspond to our derived values. 
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Figure 12. UV and far-UV observed spectra (continuous black line) of the central star of NGC 6905 and three models with T* = 150 
kK and different values of transformed radius (transformed radius and mass-loss rate are given in units of solar radii and Mq yr" 1 , 
respectively). Among the grid models, the light blue continuous line is the one that fits best the observed diagnostics, except for the two 
O V lines which are better fit by higher mass-loss rates and some weak neon features (discussed in the text). 
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Figure 13. The UV and far-UV observed spectra of the CSPN NGC 6905 (continuous black line, as in Fig |12[ | are shown along with 
three grid models with T* = 165 kK and different values of transformed radius. The discrepancy between the O V lines and all other 
diagnostics is even larger than for T* = 150 kK models (Fig. fl2[) . 
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Figure 14. The two strongest P-Cygni line profiles present in the observed spectra of NGC 6905 (continuous black line), the O VI A A 
1031.9, 1037.6 A and the C IV AA 1548.2, 1550.8 A, are compared with models with T* = 150 kK, M = 10" 7 M yr" 1 , v x = 2000 km 
s~ 1 and different turbulence velocities. 




Figure 15. The observed (continuous black line) Ly a interstellar absorption and N V AA 1238.8, 1242.8 A doublet compared with 
models with T» = 150 kK, M = 10 — 7 M.q yr — 1 , Voo = 2000 km s — 1 and different values of nitrogen abundance. 
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Figure 16. Comparison between NGC 6905 UV and far-UV observed spectra (continuous black line) and models with T* = 150 kK, 
M = 10~ 7 Mg yr — 1 , Voo = 2000 km s -1 and different values of neon abundance. The bottom panel shows the oxygen ions fractions for 
the different neon abundances. 
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Figure 17. Observed spectra of the central star of NGC 6905 (continuous black line) and models with T» = 150 kK, M = 10 7 Mq 
yr — 1 , i>oo = 2000 km s _1 and different values of oxygen abundance. The oxygen mass fraction adopted for the grid models is Xq = 0.08. 




1300 1400 1500 1600 1700 2400 2600 2800 3000 
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Figure 18. The effect on the synthetic spectra of including new ions in the calculation is seen almost exclusively in the profiles of the 
O V lines. The observed spectra are shown as a continuous black line. The pink/dashed and green/continuous lines, show models with 
and without the inclusion of new ions, respectively. The models have T, = 150 kK, M = 10~ 7 M yr _1 , v x = 2000 km s _1 . 
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Figure 19. Ar VII A 1063.55 A in models with different argon abundances. Both models shown have T, = 150 kK, M = 10 7 Mq 
yr — 1 , Voo = 2000 km s —1 . The observed spectrum of NGC 6905 is shown as a continuous black line. 
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Figure 20. The observed spectrum of the central star of NGC 6905 (continuous black line) is compared with synthetic spectra for 
different iron abundances. The models have T t = 150 kK, M = 10~ 7 M yr" 1 , v x = 2000 km s -1 . The numerous narrow absorptions 
are from interstellar H2. 




Figure 21. The Lyman-a interstellar absorption in the spectrum of NGC 6905 (continuous black line) is compared with our best-fitting 
model shown with the effect of different values of the neutral hydrogen column density. 




Figure 22. The light-blue dashed line is our best-fitting model. In pink, we apply to this model, the effect of absorption by interstellar 
molecular hydrogen with a column density of logN(H2) = 19-5 (where TV is given in units of cm -2 ). NGC 6905's far-UV observed 
spectrum is shown as a continuous black line. 
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Figure 23. NGC 6905's spectra (continuous black line) along with our best-fitting model reddened with different values of colour excess, 
assuming MW-type extinction with Ry = 3.1. 
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Figure 24. HR diagram showing the evolutionary tracks of |Miller Bertolami and Althaus] (2006 ) and the positions of our best-fitting 
grid model to NGC 6905 (which would imply a distance of 2.3 kpc), as a filled circle, and of the model rcscalcd assuming a distance of 
1.7 kpc, as an open circle. 



30 Keller et al. 



Table 1. Track A: stellar parameters of the grid models for CSPNe with 0.5Mq. 



Model 



T, [kK] log(g) R» [R ] L» [1O 3 L ] log(M) [M yr" 



Rt [R ] 



A50.M73.V500 
A50.M73.V1000 

A50.M67.V500 
A50.M67.V1000 

A50.M63.V500 
A50.M63.V1000 



500 
1000 



500 
1000 



500 
1000 



40.35 
64.04 



16.02 
25.42 



8.70 
13.80 



A65.M73.V500 
A65.M73.V1000 

A65.M67.V500 
A65.M67.V1000 

A65.M63.V500 
A65.M63.V1000 



500 
1000 



500 
1000 



500 
1000 



25.55 
40.56 



10.14 
16.10 



5.51 
8.74 



A80.M73.V500 
A80.M73.V1000 

A80.M67.V500 
A80.M67.V1000 

A80.M63.V500 
A80.M63.V1000 



500 
1000 



500 
1000 



500 
1000 



14.36 
22.80 



5.70 
9.05 



3.10 
4.91 



A100. 
A100 
A100 
A100 
A100 
A10O 
A100 
A100 
A100 



M73.V1500 
M73.V2000 
M73.V2500 
M70.V1500 
M70.V2000 
M70.V2500 
M67.V1500 
M67.V2000 
M67.V2500 



1500 
2000 
2500 



1500 
2000 
2500 



1500 
2000 
2500 



13.30 
16.12 
18.70 



8.38 
10.15 
11.78 



5.28 
6.40 
7.42 



A125.M73.V1500 125 6.3 0.08 1.5 -7.3 1500 9.43 

A125.M73.V2000 2000 11.43 

A125.M73.V2500 2500 13.26 

A125.M70.V1500 ^fji 1500 5T9T" 

A125.M70.V2000 2000 7.20 

A125.M70.V2500 2500 8.35 

A125.M67.V1500 ^T7 1500 3.74 

A125.M67.V2000 2000 4.54 

A125.M67.V2500 2500 5.26 
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Table 2. Track B stellar parameters of the grid models for CSPNe with O.GMq. 
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Model 



T, [kK] log(g) R.» [R ] L« [1O 3 L ] log(M) [M yr" 



Rt [R ] 



B50.M70.V500 
B50.M70.V1000 a 

B50.M63.V500 
B50.M63.V1000 a 

B5O.M60.V500 
B50.M60.V1000 a 



500 
1000 



500 
1000 



500 
1000 



43.80 
69.53 



14.98 
23.78 



9.45 
14.99 



B65.M70.V500 
B65.M70.V1000 a 

B65.M63.V500 
B65.M63.V1000 a 

B65.M60.V500 
B65.M60.V1000 a 



500 
1000 



500 
1000 



500 
1000 



27.86 
44.23 



9.53 
15.13 



6.01 
9.53 



B80.M70.V500 
B80.M70.V1000 a 

B80.M63.V500 
B80.M63.V1000 a 

B80.M60.V500 
B80.M60.V1000 a 



500 
1000 



500 
1000 



500 
1000 



17.59 
27.92 



6.02 
9.55 



3.79 
6.02 



B100.M70 
B100.M70 
B100.M70. 
B100.M67 
B100.M67 
B100.M67. 
B100.M65 
B100.M65 
B100.M65. 



.V1500 
.V2000 
V2500 8 
.V1500 
V2000 
V2500 8 
V1500 
.V2000 
V2500 8 



1500 
2000 
2500 



1500 
2000 
2500 



1500 
2000 
2500 



20.54 
24.88 
28.87 



12.94 
15.67 
18.19 



9.87 
11.96 
13.88 



B125.M70 
B125.M70 
B125.M70. 
B125.M67 
B125.M67 
B125.M67. 
B125.M65 
B125.M65 
B125.M65. 



.V1500 
V2000 
V2500 8 
V1500 
V2000 
V2500 8 
.V1500 
.V2000 
V2500 8 



1500 
2000 
2500 



1500 
2000 
2500 



1500 
2000 
2500 



12.95 
15.68 
18.20 



8.16 



11.46 



6.22 
7.54 
8.75 



B150.M70.V1500 
B150.M70.V2000 
B150.M70.V2500 a 
B150.M67.V1500 
B150.M67.V2000 
B150.M67.V2500 a 
B150.M65.V1500 
B150.M65.V2000 
B150.M65.V2500 a 



1500 
2000 
2500 



1500 
2000 
2500 



1500 
2000 
2500 



8.85 
10.72 
12.44 



5.58 



4.26 
5.16 
5.98 



B165.M70 
B165.M70 
B165.M70 
B165.M70 
B165.M67 
B165.M67 
B165.M67. 
B165.M67 
B165.M65 
B165.M65 
B165.M6E 
B165.M65 



.V1500 
.V2000 
V2500 8 
.V3000 
.V1500 
.V2000 
V2500 8 
.V3000 
V1500 
.V2000 
V2500 8 
V3000 



1500 
2000 
2500 
3000 



1500 
2000 
2500 
3000 



1500 
2000 
2500 
3000 



6.67 



9.37 
10.58 



4.20 
5.09 
5.90 
6.67 



3.21 
3.88 
4.51 
5.09 



B200.M73 
B200.M73 
B200.M73 
B200.M73 
B200.M70 
B200.M70 
B200.M70 
B200.M70 
B200.M67 
B200.M67 
B200.M67. 
B200.M67 
B200.M65 
B200.M65 
B200.M65. 
B200.M65 



V1500 
.V2000 
V2500 
.V3000 
.V1500 
.V2000 
V2500 a 
.V3000 
V1500 
.V2000 
V2500 a 
.V3000 
.V1500 
.V2000 
V2500 a 
V3000 



1500 
2000 
2500 
3000 



1500 
2000 
2500 
3000 



1500 
2000 



3000 



1500 
2000 
2500 
3000 



4.60 
5.57 
6.47 
7.31 



2.90 
3.51 
4.08 
4.6 



1.83 
2.21 



2.90 



1.39 
1.69 
1.96 
2.21 
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Table 3. Track C: 


stellar parameters of the grid models for CSPNe with 0.9Mq . 






Model 


1* [kKJ log(g) R* [R©J L* [10 Lq\ log(M) [M©yr J 


-It 

3 [km s ] 


Rt [R©] 


C100.M67.V1500 


100 4.8 0.62 33.5 -6.7 


1500 


27.59 


C100.M67.V2000 




2000 


33.42 


C100.M67.V2500 




2500 


38.79 


C100.M65.V1500 


-6.5 


1500 


21.06 


C100.M65.V2000 




2000 


25.51 


C100.M65.V2500 




2500 


29.60 


C100.M64.V1500 


-6.4 


1500 


17.38 


C 100. M64.V2 00(1 




2000 


21.06 


C100.M64. V2500 




2500 


24.44 


C125.M67.V1500 


125 5.3 0.35 25.8 -6.7 


1500 


15.47 


C125.M67.V2000 




2000 


18.75 


C125.M67.V2500 




2500 


21.75 


C125.M65.V1500 


-6.5 


1500 


11.81 


C125.M65.V2000 




2000 


14.31 


C125.M65.V2500 




2500 


16.60 


C125.M64.V1500 


-6.4 


1500 


9.75 


C125.M64.V2000 




2000 


11.81 


C125.M64.V2500 




2500 


13.70 


C165.M67.V1500 


165 5.7 0.22 31.2 -6.7 


1500 


9.75 


C165.M67.V2000 




2000 


11.82 


C165.M67.V2500 




2500 


13.71 


C165.M67.V3000 




3000 


15.48 


C165.M65.V1500 


-6.5 


1500 


7.44 


C165.M65.V2000 




2000 


9.03 


C165.M65.V2500 




2500 


10.46 


C165.M65.V3000 




3000 


11.82 


C165.M64.V1500 


-6.4 


1500 


6.14 


C165.M64.V2000 




2000 


7.44 


C165.M64.V2500 




2500 


8.64 


C165.M64.V3000 




3000 


9.75 


C200.M67.V1500 


200 6.0 0.15 33.8 -6.7 


1500 


6.90 


C200.M67.V2000 




2000 


8.36 


C200.M67.V2500 




2500 


9.71 


C200.M67.V3000 




3000 


10.96 


C200.M65.V1500 


-6.5 


1500 


5.27 


C200.M65.V2000 




2000 


6.38 


C200.M65.V2500 




2500 


7.41 


C200.M65.V3000 




3000 


8.36 


C200.M64.V1500 


-6.4 


1500 


4.35 


C200.M64.V2000 




2000 


5.27 


C200.M64.V2500 




2500 


6.11 


C200.M64.V3000 




3000 


6.90 
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Table 4. Species considered in the calculation of the grid models. 



clement 


ions 


ne 


T TT TTT 

lj 11, 111 




TTa TTTa T\/ V 
11 , 111 ,1V, V 


N 


TTa jTTa j\r& sj \/T 
11 > 111 ,1V , V , VI 


O 


IP, IIP, IV a , V, VI, VII 


Ne 


IP, IIP, IV a , V, VI, VII, VIII, IX 


Al 


IIP, IV a , v a 


Si 


IIP, IV, V a , VP 


P 


IV a , V, VI 


S 


IIP, IV a , V a , VI, VII 


Fc 


IV a ya VI a VII) VIII; IX; X, XI 



a These ionic species are not included in the calculation of all models. They were introduced where needed based on the analysis of the 
ionization fractions. 
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Table 5. Spectra of NGC 6905's central star used in the analysis. 



Instrument 


Data Set 


Date 


Resolution [A] 


Aperture [arcsec] 


Range [A] 


FUSE 


A1490202000 


2000 Aug 11 


~0.06 


30x30 


905-1187 


STIS + G140L 


O52R01020 


1999 Jun 29 


~1.20 


52x0.5 


1150-1736 


STIS + G230L 


O52R01010 


1999 Jun 29 


~3.15 


52x0.5 


1570-3180 



